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The Hall and Seebeck coefficients, electrical resistivity and thermal conductivity of polycrystalline 
NiSiaP4 were characterized from 2 to 775 K. Undoped NiSiaP4 behaves like a narrow gap semicon- 
ductor, with activated electrical resistivity p below room temperature and a large Seebeck coefficient 
of ~ 400/iV/K at 300 K. Attempts to substitute boron for silicon resulted in the production of ex- 
trinsic holes, yielding moderately-doped semiconductor behavior with p increasing with increasing 
temperature above ~ 150 K. Hall carrier densities are limited to approximately 5xl0 19 cm -3 at 
200 K, which would suggest the solubility limit of boron is reached if boron is indeed incorporated 
into the lattice. These extrinsic samples have a Hall mobility of ~12cm 2 /V/s at 300 K, and a 
parabolic band equivalent effective mass of ~3.5 times the free electron mass. At 700 K, the ther- 
moelectric figure of merit zT reaches ~0.1. Further improvements in thermoelectric performance 
would require reaching higher carrier densities, as well as a mechanism to further reduce the lattice 
thermal conductivity, which is ~ 5 W/m/K at 700 K. Alloying in Ge results in a slight reduction of 
the thermal conductivity at low temperatures, with little influence observed at higher temperatures. 
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I. INTRODUCTION 

Identifying thermoelectric materials composed of en- 
vironmentally friendly and/or plentiful elements is one 
current thrust of the thermoelectric community. Recent 
successes in this area includ e se veral Zintl phases, such 
as Ca 3 AlSb3 and CasAUjSbeP^l Simple binaries such as 
FeSb2 and FeSi may be promising for low-temperature 
applications^^ In addition to environmental concerns 
and abundance issues, low-density materials are also of 
interest for applications due to their potential for large 
specific powers, and the associated reduced mass require- 
ments. For instance, Mg2Si is an example of a material 
that possesses all of the above desired properties and 
good thermoelectric performance!^! The higher man- 
ganese silicides (MnSi2-a;, a;~0.25) are another group of 
low-cost, abunda nt m aterials with promising thermoelec- 
tric performanceP^ It is this combination of abundant 
and low-density materials that brought us to examine an- 
other silicide, NiSi3P4, of which little is known beyond 
the crystallographic structure. 

NiSiaP4 is isostructural with Cu3SbSe4p^ which has 
recently been shown to possess moderate thermoelectric 
performance. At the materials level, thermoelectric per- 
formance is quantified via the dimensionless figure of 
merit zT=a 2 T / pn, where a is the Seebeck coefficient, 
p the electrical resistivity, and k the thermal conductiv- 
ity. In Cu3SbSe4, zT reaches between 0.7 and 0.9 near 
^650 K, with the exact value depending on the particular 
alloying and doping employed E2H2D 

NiSi3P4 and Cu3SbSe4 are isostructural with the min- 
eral famatinite Cu3SbS4, which is a member of the stan- 
nite (Cu2FeSnS4) group. These materials possess a 
tetragonal unit cell (space group /42m, No. 121) based 
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on the zinc-blende structure with all atoms being tetrahc- 
drally coordinated (though the tetrahedral coordinations 
are slightly distorted with bond angles deviating from the 
ideal 109.5°). In Cu 3 SbSe4, Cu and Sb are tetrahedrally 
coordinated by Se, whereas in NiSi3P4, Ni and Si are 
surrounded by P. 

In this work, we characterize NiSisP4, and find it to 
be a narrow gap semiconductor where holes naturally 
dominate conduction. This material is too resistive to 
yield large thermoelectric performance, and thus doping 
studies were undertaken. At first, we attempted to sub- 
stitute Co for Ni, though this did not appear to produce 
extrinsic holes. Thus, in analogy with Si and SiGe alloys, 
we attempted to use boron as a hole dopant, and nom- 
inal compositions NiSi 3 _ x B a; P4 were found to produce 
extrinsic hole concentrations. If boron substitutes for sil- 
icon, it does so with a relatively limited solubility, as the 
maximum hole concentration observed in these samples 
is approximately 5xl0 19 cm -3 at 200 K. It is also pos- 
sible the extrinsic hole concentration is associated with 
off-stoichiometry of the NiSisP4 phase induced by the 
nominal compositions employed. With some modest suc- 
cess at doping, it was clear that the lattice thermal con- 
ductivity severely limits the thermoelectric performance. 
Thus, we attempted to alloy in 20% Ge on the Si site, 
and this resulted in only a small decrease in k (primarily 
at low temperatures). The net result is a rather low fig- 
ure of merit zT^O.l at 700 K. Further improvements in 
performance may be achieved with higher doping levels, 
as the material has a relatively high band mass, though a 
reduction in the lattice thermal conductivity is also nec- 
essary. Therefore, this material is unlikely to be pursued 
in further detail unless theoretical calculations provide 
insights into doping or suggest anomalous performance 
for other regions of the electronic structure. 
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II. METHODS 

Polycrystalline samples were prepared by first arc- 
melting high-purity Ni, Si, B, and Ge as appropriate for 
the nominal compositions listed in Table [TJ The resulting 
ingot was hand ground to a fine powder, and combined 
with red phosphorus in a silica ampoule that was subse- 
quently sealed under vacuum. The ampoule was heated 
to 450 °C, held at this temperature for 10 h, then heated 
to 900 °C and held for 100 h to allow phosphorus to re- 
act through the vapor phase. After cooling, the ampoule 
was shaken vigorously and heated to 900 °C for another 
^50-100 h to ensure full reaction of the phosphorus. Dur- 
ing the first heating, a rate of 25°C/h was employed, 
and 50°C/h was used for the second heat treatment. 
Indications of decomposition and significant volatiliza- 
tion of phosphorus were observed when sintering of a 
cold-pressed pellet at 1000 °C was attempted, and thus 
^900 °C was used as a high temperature limit. 

After heating, the powders were hand ground and 
consolidated with a hot press. Hot pressing was per- 
formed in a graphite die within a graphite furnace, us- 
ing pressures near 10,000 psi and a maximum tempera- 
ture near 900 °C±50 °C. After exposure to the maximum 
force/temperature for ^lh, the force was removed for a 
one hour stress- free anneal, followed by cooling to 675 °C 
over one hour, at which temperature the furnace was 
turned off to cool more rapidly. This procedure resulted 
in samples with densities greater than 92% of the crys- 
tallographic density. The samples are relatively stable in 
air, with only a very slight tarnish appearing after several 
months. Grinding was generally done in a helium glove 
box, though grinding in air prior to hot pressing did not 
appear to affect the transport properties provided the re- 
action was complete. The acrid smell of phosphorus was 
experienced upon cutting or grinding the final specimens. 

The polycrystalline samples were found, by powder x- 
ray diffraction, to contain mainly the desired NiSiaP4 
phase. The undoped NiSi3P4 sample contained a small 
amount (<2 wt.%) of N12SL Rietveld refinement of the 
room temperature powder diffraction data using the pub- 
lished structural model 11 yielded lattice parameters of 
a=5.155(l)A and c=10.345(l)A. Boron substitution re- 
sulted in a negligible change in lattice parameters. With 
increasing nominal concentrations of boron, the impu- 
rity phase(s) changed from N^Si to BP and NiSi, with 
up to 3% and 6.4%, respectively, observed in sample 4,B 
in Table H] As discussed below, transport measurements 
suggest that the boron solubility limit was reached by 
the nominal composition of NiSi2.9sB0.02P4- It is also 
possible that the phase width of NiSi3P4 is reached at 
this nominal composition and boron does not incorpo- 
rate into the lattice. In the sample with nominal com- 
position NiSi2Geo.6Bo.4P4, the diffraction revealed about 
87% purity of the main phase, with BP (7 wt.%) and Ge 
(6 wt.%) as the primary secondary phases. The lattice 
parameters indicate some Ge incorporated into the lat- 
tice, though, with refinement yielding a=5.169(l)A and 



TABLE I. Nominal compositions and lattice parameters of the 
N1S13P4 and doped-NiSi3P4 samples presented in this study, 
sample ID nominal composition a c 



1 NiSi 3 P 4 

2, B NiSi2.9sB0.02P4 

3, B NiSi2.90B0.10P4 

4, B NiSi2.80B0.20P4 

5, B & Ge NiSi2Ge0.eB0.40P4 



A A 

5.155(1) 10.345(1) 
5.155(1) 10.345(1) 
5.157(1) 10.347(1) 
5.154(1) 10.342(1) 
5.169(1) 10.367(1) 



c=10.367(l)A. 

Transport data suggested very high hole concentra- 
tions are necessary to achieve large thermoelectric per- 
formance, and thus we attempted to synthesize sam- 
ples with large concentrations of boron. The secondary 
phases seem to have little influence on the transport 
properties, as similar results were obtained on the dif- 
ferent samples. To allow for an appreciation of the varia- 
tions in properties across batches, we report data on sev- 
eral boron containing samples that possess similar trans- 
port properties and thus likely possess similar amounts 
of boron in the primary phase. 
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FIG. 1. (color online) (a) Electrical resistivity, (b) Seebeck 
coefficient, and (c) thermal conductivity of polycrystalline 
MS13P4, as well as B- and Ge-doped compositions as listed 
in Table HI 

Seebeck coefficient a, electrical resistivity p, and Hall 
coefficient Rh measurements were utilized to investi- 
gate the electrical properties; the Hall carrier density 
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is nH—l/Rne. Below 300 K, these measurements were 
performed in a Quantum Design Physical Property Mea- 
surement System (PPMS). The thermoelectric measure- 
ments, including the thermal conductivity n, were per- 
formed using the PPMS Thermal Transport Option 
(TTO), and the Hall effect was measured using the re- 
sistivity option on the PPMS. For TTO measurements, 
gold-plated copper leads were attached to the sample 
using silver epoxy (H20E Epo-Tek), and silver paste 
(DuPont 4929N) was utilized to obtain low resistance 
contacts for Hall effect measurements performed using 
0.025 mm Pt wires. Hall coefficients were obtained from 
a fit of the Hall resistance versus magnetic field, with 
maximum fields of ±6 Tesla employed. 

High temperature thermoelectric performance was as- 
sessed using an ULVAC ZEM-3 M8 to measure a and p. 
The thermal conductivity k (above room temperature) 
was calculated via k = Dt x Cp x d where d is the den- 
sity. The specific heat capacity C p was approximated 
using the Dulong Petit limit of C v = 3R, calculated ac- 
cording to nominal composition. The differences between 
C v and Cp are assumed to be negligible at the temper- 
atures of interest, and this estimation likely results in a 
slight overestimation of zT at the highest temperatures. 
The thermal diffusivity Dt was obtained on ~10mm di- 
ameter discs between 1.5 and 2 mm thick, in an Anter 
FL-5000 with a graphite furnace; data analysis followed 
ASTM 1461 for flash diffusivity. No sign of oxidation oc- 
curred during the high temperature measurements, and 
data collected on cooling were consistent with those col- 
lected on heating. 



III. RESULTS AND DISCUSSION 

NiSi3P4 was found to behave similar to a narrow gap 
semiconductor with a low concentration of charge carri- 
ers. As shown in Figure [T] activated electrical conduc- 
tion is observed for MS13P4 below room temperature. 
Between 150 and 250 K, a simple activation energy was 
obtained from fitting to p= J ooExp[£ , Q /kT], resulting in 
i? a =0.11eV, which is equivalent to a fundamental band 
gap of 0.22 eV. However, the fit to this simple activated 
behavior becomes poor when wider temperature ranges 
are considered. 

Consistent with the semiconducting-like behavior, a 
large Seebeck coefficient is observed in NiSiaP4, as shown 
in Figure [TJd. Also consistent with this behavior, the Hall 
carrier density uh increases with increasing temperature 
for NiSisP4, as shown in Figure [2^,. We emphasize that 
this is a single band n# , and may not reflect the true car- 
rier density if electrons are also contributing to electrical 
transport. 

In order to probe the potential for thermoelectric per- 
formance, we attempted to hole-dope NiSisP4. Wc first 
tried the substitution of Co for Ni, but this sample had 
a higher resistivity than undoped NiSi3P4 at 300 K, and 
also displayed activated conduction below room temper- 



ature. We then found that boron successfully increased 
the hole concentration in samples with nominal composi- 
tions of NiSi3_ :r B a ,P4. However, synthesis with £=0.02, 
0.10, 0.20, 0.40 all yielded similar transport properties, 
as observed in Figure [T] One explanation for this trend 
is that the solubility limit of boron has been reached. 

Hole-doped samples of NiSi3_ a; B a ;P4 behave metallic 
above ~ 150 K with p increasing with increasing T. Sim- 
ilarly, the Seebeck coefficients of these samples increase 
with increasing T, consistent with extrinsically doped 
samples. The Seebeck coefficients remain rather large 
in these doped samples, with a~250/iV/K obtained at 
300K. 

The Hall carrier concentration and mobility pn of 
boron- and boron/germanium-doped MS13P4 are shown 
in Figure [2] As observed, the values of nn are relatively 
temperature independent, and fairly similar for the dif- 
ferent nominal boron concentrations. This is consistent 
with similar p and a in the various samples. This re- 
sult suggests that, if boron substitutes for silicon, the 
solubility limit of boron has been reached. If lattice- 
incorporated boron atoms are fully ionized and con- 
tribute one hole to the system, then a solubility limit of 
^0.2%B per silicon is obtained from the Hall carrier con- 
centration of approximately 5x 10 19 holes/cm 3 at 200 K. 
This seems reasonable considering the solubility limits 
of boron in silicon, which are ~ 2xl0 19 B/cm 3 at 700°C 
and - lxl0 20 B/cm 3 at 900° CP The small increase of 
tih above ~200K may be the result of some carriers be- 
ing activated from other impurities or across the band 
gap. 

The hole mobility of boron-doped NiSiaP4 is found 
to be approximately 12cm 2 /V/s at 300 K. The mobility 
increases with decreasing temperature, consistent with 
a carrier mean free path limited by phonon scattering. 
However, this behavior only occurs down to ^125 K, 
where the mobility reaches a maximum of approximately 
18cm 2 /V/s. While pn decreases upon cooling below 
125 K, tih remains relatively constant, which may sug- 
gest that additional carrier scattering mechanisms are 
important at low T . Ionized impurity scattering gener- 
ally leads to a mobility that increases as T 3 / 2 , though 
the low-T rise in pn occurs more slowly than T 3 / 2 in 
boron-doped NiSisP4. It is also possible that Anderson 
localization occurs due to disorder in the lattice, causing 
a mobility gap that leads to an increase in p at low- 
T. Also, we note that the high-temperature decay oc- 
curs more slowly than the T -3 / 2 behavior expected for 
a lightly doped semiconductor (or even the T _1 behav- 
ior expected for a degenerate semiconductor). Therefore, 
it is possible impurity band conduction is important in 
these materials or grain boundary scattering influences 
transport. 

The hole mobility of boron-doped NiSi3P4 is less than 
that of isostructural Cu3SbSe4. The mobility of undoped 
Cu 3 SbSe 4 is ~41cm 2 /V/s at 300 K, with a hole con- 
centration of 7.5x 10 18 cm~ 3 P The mobility of undoped 
NiSisP4 is not shown in Figure [2] because the single band 
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FIG. 2. (color online) (a) Hall carrier concentration and (b) 
Hall mobility of polycrystalline NiSisP4, as well as B- and 
Ge-doped compositions as listed in Table [I] 



model likely fails for the undoped composition. At a 
more comparable hole concentration of 6.3xl0 19 cm -3 
in Cu3Sbo.975Sno.o25Se4 the mobility is approximately 
35cm 2 /V/s, and it is 31cm 2 /V/s in Cu 3 Sbo.95Sn .o5Se4 
at 300 KP^ Given that the carrier mobility in doped 
NiSiaP4 is dominated by phonon scattering at room tem- 
perature, it seems unlikely that the boron dopant is lead- 
ing to this large difference in carrier mobility. Therefore, 
the hole mobility in NiSisP4 appears to be inherently 
lower than in the copper-based analogues. 

To better understand the differences between trans- 
port in NiSi3P4 and Cu3SbSe4, we consider a parabolic 
band analysis of the Hall and Seebeck coefficients (Equa- 
tions [lJQ . Analysis of the Hall and Seebeck coefficient 
data allows the carrier mass m* to be obtained. When 
phonon scattering limits the relaxation time of charge 
carriers, the mobility is expected to scale as /^cxto*~ 5 / 2 . 
Therefore, even small changes in the band mass can have 
significant influence on carrier transport. 

Within a single parabolic band model, the Seebeck co- 
efficient is a function of the reduced electrochemical po- 
tential r\ and scattering parameter A: 



k f {2 + \)F x+1 {r 1 ) 
(1 + A)F A (77) 



n 



(i) 



The scattering parameter A is used to incorporate the 
energy dependence of the carrier relaxation time, r = 
Toe A ~ ' 5 . The value of A influences a, but To does not. 
In this analysis, we assume acoustic phonon scattering 
limits r, corresponding to A=0. 

The Hall carrier density, nn=\j i?#e where Rh is the 
Hall coefficient, is given by 
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FIG. 3. (color online) High temperature (a) electrical resis- 
tivity, (b) Seebeck coefficient, and (c) thermal conductivity of 
polycrystalline NiSisP4, as well as B- and Ge-doped compo- 
sitions as listed in Table [U 
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Analysis of nn allows the effective mass to* to be deter- 
mined for a given 77 obtained from analysis of the Seebeck 
coefficient. Carrier scattering influences nn through the 
Hall factor r#, which is a function of r\ and A similar to 
the Seebeck coefficient: 



3 ^ (l/2 + 2A)F 2A _ 1/2 (r ? ) 

r "=2 F ^ (l + A) 2 F 2 (r,) ■ (3) 

For simplification purposes, these expressions utilize the 
Fermi integrals Fj{rf) defined by 



1 + Ex P [£ - rj\ ' 



(4) 



Employing the above equations, we find that the ef- 
fective mass is approximately 3.5±0.5m e in boron-doped 
NiSisP4 at room temperature; m e is the free electron 
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mass. At 200 K, where nn is still relatively independent 
of temperature, this process yields 2.9-3.7m e for the var- 
ious samples of boron-doped NiSiaP4. From literature 
values for nn and a in Cu3SbSe4 (corresponding to the p, 
values noted above) p21 we obtain 0.45m e < m* < 1.3m e 
depending on the particular sample. Therefore, it is clear 
that the band mass is much larger in NiSi3P4 than in 
Cu3SbSe4, and this likely contributes to the lower mo- 
bility. This large band mass results in a relatively large 
Seebeck coefficient in NiSi3P4, but the impact on carrier 
mobility is detrimental to the thermoelectric performance 
ofNiSi 3 P 4 . 

The thermal conductivity of NiSi3P4 is approximately 
lOW/m/K at 300K. Boron doping reduces the low- 
temperature value of K slightly. In an attempt to fur- 
ther reduce k, we substituted Ge for Si, in addition to 
the boron doping. This results in some additional sup- 
pression of k at low-T, though little change is observed 
at higher temperatures. Interestingly, k does not display 
much temperature dependence below room temperature, 
even in undoped NiSi3P4. This suggests there may be 
some disorder in the system, or perhaps grain bound- 
aries or secondary phases are reducing n at low-T. We 
note that due to the large values of p, the lattice ther- 
mal conductivity kl dominates k for all samples below 
~600K. 

The high temperature transport properties are plotted 
in Figure [3] Interestingly, the decrease in p with increas- 
ing T for undoped NiSiaP4 turns into a nearly temper- 
ature independent p above ~400 K, and the correspond- 
ing Seebeck coefficient is large and relatively tempera- 
ture independent. The boron-doped samples display typ- 
ical doped semiconductor behavior at high temperatures, 
with the resistivity and Seebeck coefficients increasing 
with increasing temperature until the effects of minor- 
ity carrier activation are observed near 650 K. The ther- 
mal conductivity decreases with increasing T in all sam- 
ples, consistent with the dominant role of phonon-phonon 
scattering in these crystalline materials. As noted above, 
the substitution of 20% Ge (nominal composition) did 
not result in a significant reduction in n above 300 K. For 
comparison, the thermal conductivity of Cu3SbS4 reaches 
approximately 1.5W/m/K at 600 KP 

The maximum in a(T), which is relatively unpro- 
nounced in these materials, can be utilized to estimate 
the thermal band gap by E g =2a max T max W^ This anal- 
ysis suggests that the band gap is between 0.43 and 
0.50 eV, which is about two times larger than that ob- 
tained from the analysis of p(T) at lower temperatures. 
It is possible that one of these estimates of E g is influ- 
enced by the activation of carriers from an impurity band. 
Without high temperature Hall effect data the source for 
this discrepancy is difficult to isolate. 

The thermoelectric power factor, a 2 T/p is plotted 
in Figure [4] along with the figure of merit zT. Rel- 
atively low thermoelectric performance is observed in 
these materials. This is due in part to the relatively high 
thermal conductivity for the corresponding mobility val- 
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FIG. 4. (color online) (a) The thermoelectric power factor 
times temperature and (b) dimensionless figure of merit for 
polycrystalline NiSi3P4, as well as B- and Ge-doped compo- 
sitions as listed in Table |U 



ues. Boron-doping greatly improves both the power fac- 
tor and zT , and additional improvements would likely 
come with higher doping levels. This is particularly 
the case for p-type NiSisP4, because it is a large band 
mass material and the optimum carrier density increases 
with increasing band mass and increasing lattice thermal 
conductivity!^^ As noted, however, boron substitution 
only permits modest carrier densities to be reached, and 
using a transition metal dopant was not successful in in- 
creasing the hole concentration. It is possible that al- 
tering the Si/P ratio would permit higher doping levels 
to be achieved, and may also allow n-type behavior to 
be probed. However, it seems unlikely that this material 
will achieve the zT values observed in isostructural cop- 
per compounds, which have higher mobility due to lower 
band masses, as well as lower lattice thermal conductiv- 
ities. 



A. Summary 

Polycrystalline NiSiaP4 was found to be a narrow gap 
semiconductor, with an activation energy for the electri- 
cal resistivity on the order of lOOmeV below room tem- 
perature. Boron-doping is successful in producing ex- 
trinsic holes, though the corresponding hole concentra- 
tion is limited to approximately 5xl0 19 cm~ 3 , suggest- 
ing a possible solubility limit of boron in NiSi3P4. The 
transport properties of the boron-doped samples behave 
as expected for moderately doped semiconductors, and 
parabolic band analysis suggests a fairly high band mass 
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of ^3m e . The thermoelectric figure of merit only reaches 
~0.1 at 700 K, due in large part to a lower than desired 
hole mobility and carrier concentration, and a moder- 
ately large thermal conductivity. Further improvements 
in thermoelectric performance may be possible with a 
different doping scheme or through additional phonon 
scattering, but are unlikely to yield large thermoelectric 
efficiencies. 
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